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Abstract: Perturbative gravity about global de Sitter space is subject to linearization- 
stability constraints. Such constraints imply that quantum states of matter fields couple 
consistently to gravity only if the matter state has vanishing de Sitter charges; i.e., only 
if the state is invariant under the symmetries of de Sitter space. As noted by Higuchi, 
the usual Fock spaces for matter fields contain no de Sitter-invariant states except the 
vacuum, though a new Hilbert space of de Sitter invariant states can be constructed via 
so-called group- averaging techniques. We study this construction for free scalar fields of 
arbitrary positive mass in any dimension, and for linear vector and tensor gauge fields in 
any dimension. Our main result is to show in each case that group averaging converges for 
states containing a sufficient number of particles. We consider general A^-particle states 
with smooth wavefunctions, though we obtain somewhat stronger results when the wave- 
functions are finite linear combinations of de Sitter harmonics. Along the way we obtain 
explicit expressions for general boost matrix elements in a familiar basis. 
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1. Introduction 

It has been known since the 1970s [1, 2, 3, 4, 5, 6] that perturbative gravity about a back- 
ground having both compact Cauchy surfaces and Killing symmetries is subject to so-called 
linearization-stability constraints; i.e., that a linearized perturbation can be extended to a 
full solution only of it satisfies certain constraints. These constraints may be summarized 
by noting that linearized gravity on such backgrounds has conserved charges corresponding 
to the Killing symmetries. In the classical context, the constraints state that a solution 
of linearized gravity can be extended to a solution of the full interacting theory only if 
these conserved charges vanish. In the quantum context, the vanishing of such charges 
requires the Hnearized state to be invariant under the symmetries of the background [7] . 
When matter fields are present, it is the state which describes both matter and linearized 
gravitational waves propagating on the fixed background which much be invariant. The 
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physical root of this phenomenon is that KiUing symmetries of the background are gauge 
symmetries when Cauchy surfaces are compact. 

Perturbative quantum gravity about de Sitter space presents a particularly interesting 
case of this phenomenon. Perturbative states must be invariant under the de Sitter group 
[8, 9, 10] but, because this group is non-compact, the usual Fock spaces contain no de 
Sitter invariant states except for a possible vacuum state [8]; i.e., they contain far too few 
states for a useful theory of perturbative gravity. 

However, it was shown in [9] that one may use a standard Fock space of non-invariant 
states (which we call the 'auxiliary' space Tia.ux) to construct a new 'physical' Hilbert space 
"^phys of de Sitter invariant states via techniques known as group averaging. The issues 
involved are very similar to those that arise in Dirac quantization of constrained systems 
[11]; indeed, group-averaging was independently introduced in that context [12, 13]. The 
basic idea is to note that linear combinations of auxiliary states of the form 



Jg&G 

are formally de Sitter invariant. Here U{g) gives the unitary representation of G on Waux 
and dg is the Haar measure on the de Sitter group G. Since the de Sitter group is unimodu- 
lar, dg is invariant under both right- and left-translations. For compact groups the integral 
(1.1) converges and simply projects the state onto the trivial representation of G. For 
non-compact G the state (1.1) is not normalizeable in TCaux, but can nevertheless often be 
understood [14, 15] as a well-defined "generalized state" in a sense similar to that used 
in quantum mechanics for the non-normalizeable eigenstates of operators with continuous 
spectrum (e.g., plane waves in infinite space). The key questions are whether the so-called 
group averaged inner product 



converges on a sufficiently large set of states, and whether this inner product is positive 
definite. If so, it make be taken^ as the physical inner product which defines TCphys- 

The group- averaged inner product (1.2) and the Hilbert space Hphys were studied in 
detail for the special cases of free 2d massless scalar fields (with the zero-mode removed) 
and 4d linearized gravitons in [9] and for 4d conformally coupled scalars in [19]. For these 
cases it was shown that (1.2) converges when the states 1^/^1,2) contain sufficiently many 
particles and that it is positive definite. For these special cases it was possible to find a 
fairly explicit orthonormal basis for TCphys- (See also [20] for a similar treatment including 
the zero-mode for the 2d massless scalar.) 

However, the inner product (1.2) is more difficult to study for general free fields. 
Such an analysis was begun in [19], which showed that group averaging converges for 1-1-1- 
dimensional principal series massive scalar fields in certain states with > 3 particles. Our 

^The point of this construction is that, since aU gauge-invariant observables commute with U{g), they 
have a natural action on Hphys- For comments on the extraction of physics from such observables see 




(1.1) 





[7, 13, 16, 17, 18]. 
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purpose here is to generalize the known results to arbitrary dimensions, arbitrary states 
with smooth n-particle wavefunctions, and more general fields. In particular, we consider 
arbitrary scalar masses > as well as both vector gauge fields and linearized gravitons. 
Our main result will be that (1.2) converges whenever the auxiliary states IV'1,2) contain 
suflttciently many particles for either i) general smooth A^-particle wavefunctions or ii) the 
special case of wavefunctions given by finite linear combinations of de Sitter Harmonics. For 
linear vector gauge fields, linearized gravitons, and scalar fields with masses corresponding 
to principal series representations of the dSd+i de Sitter group, d + 1 particles suffice for 
general smooth wavefunctions. For scalar fields of smaller mass, we require more particles 
in a way that will be explained after establishing our conventions below. For most cases, 
we are able to show convergence with even fewer particles for the special case (ii). As 
discussed in section 4, we suspect that these smaller particle numbers may also suffice for 
the more general case (i). 

While our analysis does not guarantee that (1.2) is positive definite, our results do 
show that (1.2) is positive when jipi) = lip2) is a tensor product of our de Sitter harmonics. 
In addition, a theorem of [21] states that when (1.2) converges it gives the unique inner 
product on states of the form (1.1) consistent with the *-algebra of bounded gauge-invariant 
observables on Waux- Thus on physical grounds one expects that (1.2) is positive definite^. 

Our discussion is organized as follows. Section 2 examines group averaging for massive 
scalar fields and obtains explicit expressions for the matrix elements of U{g) in a familiar 
basis, though certain technical details are relegated to the appendices. Such expressions 
allow us to study the convergence of (1.2) for scalar fields, and to establish that group 
averaging converges under the conditions stated above. We then perform a similar study 
for vector and tensor gauge fields in section 3. We find that many results from the scalar 
case can be applied directly to the vector and tensor cases. We conclude with a brief 
discussion of open issues and states with small particle numbers in section 4. 

2. Massive scalar fields 

We wish to analyze the convergence of group averaging for massive scalar fields in D = d+1 
dimensional de Sitter space, generalizing and extending certain results of [19] for the case 
D = 2. Our main task will be to construct an explicit expression for the matrix elements 
{tpilB{X)l'ip2) for states IV'1,2) in the standard de Sitter Fock space of a free massive scalar 
field, from which the convergence properties of (1.2) can be easily evaluated. Here B{\) is 
a finite boost of rapidity A. As discussed in section 2.3 below, such matrix elements control 
the convergence of (1.2). 

In an A^-particle state, boost matrix elements reduce to (symmetrized) products of 
matrix elements for one-particle states. We will therefore focus on computing ('i/'i|i?(A)|'02) 
for one-particle states below. These are just matrix elements for irreducible representations 
of the de Sitter group. There is a long history of such computations in the literature, 
beginning with Bargmann's work [25] on S0(2,l) and continued in [26, 27] for S0(3,l); see 

^See [22, 23, 24] for further comments on the positivity of (1.2) in other settings. 
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also [28, 29]. These works compute matrix elements in the basis we use below. Using the 
rather different Gelfand-Tsetlin basis [30], a general formula for all S0(-D,1) was found in 
[31] and was evaluated explicitly in [32] in terms of hypergeometric functions. Appendix B 
is devoted to providing a more direct computation of the general S0(Z),1) matrix elements 
for scalar representations in a basis more familiar to physicists. For our purposes the results 
of appendix B are somewhat simpler than those of [32] and give us better control over the 
asymptotics at large g. 

We establish notation and review scalar fields on de Sitter in section 2.1. Section 2.2 
computes the desired matrix elements of B(X) using results from the appendices. The 
convergence of (1.2) is then analyzed in section 2.3. 

2.1 Scalar Fields in de Sitter 

Recall that the de Sitter metric may be conveniently described by embedding dSo in a 
Z) + 1-dimensional Minkowski space M/j+i. Using Cartesian coordinates X"^ in Mjj^i, one 
may define de Sitter spacetime as the hyperboloid of constant radius £: 

e = t^abX^X'' = -{Xy + (X^ + ■■■ + {X'^f , (2.1) 

where r]AB is the D + 1 Minkowski metric. The induced metric can be written 

ds^ = e [-dt^ + (cosh^ t)d^fJ\ 

= f [-dt^ + (cosh^ t) {dx^ + sin^ x dl^Li)] , (2.2) 

where t = si'Dh~^{X^ /£), d = D — 1, and is the metric of the d-sphere S"^. 
We consider scalar fields satisfying the Klein-Gordon equation 

'a{a + dy 



□(^ = (M^ + ^R) if 



£2 



^ , (2.3) 



where is the mass, R = D{D — 1)/^^ is the scalar curvature of de Sitter, ^ is a non- 
minimal coupling constant, and □ := g^'^V^'Vu is the wave operator on de Sitter. We define 
the complex parameter a through the last equality. Below, we focus on the case a ^ 0, —d, 
though this case can be included in the analysis below if the zero-mode is removed. (See 
[20] for a discussion of group-averaging the a = zero-mode for D = 2.) 

One may quantize the scalar field (pa with Casimir —a{a + d) by expanding it in a sum 
over modes 

¥'a(x) = ajip^j{x) + aijp^jix) , (2.4) 
/ 

where the ip^j form an irreducible representation of the de Sitter group. The vector j = 
jd) • • • ) Ji contains d angular momentum quantum numbers which satisfy 

id >id-i > ••• > |ji| ;. (2.5) 

Below, it will be convenient to denote the total angular momentum by j := jd and the 
"second" quantum number by := jd~i- See [33, 34, 35, 36] for classic references on 
quantizing such fields. 



-4- 



We take our modes if ^ to he normalized with respect to the Klein-Gordon inner 
product 

so that the creation and annihilation operators satisfy the usual algebra: 



6^^ 

jm 



at, al 



. (2.7) 



They can also be shown to satisfy a de Sitter invariant positive frequency condition (see 
e.g. [37]). 

The auxiliary Hilbert space Waux is then built from the vacuum []0) satisfying 

a.lO) = Vj. (2.8) 

The discussion above is explicitly de Sitter-invariant, so |0) is the familiar de Sitter- 
invariant (Bunch-Davies) vacuum. The full auxiliary space Ti^ux is spanned by states 
of the form lip) = at ... at |0), whose inner products are readily computed from (2.6) and 
(2.7). 

The one-particle states of Waux form a scalar representation of S0q{D,1). Recall 
[38, 25, 39, 40, 41] that scalar representations T^'^ of SOo{D, 1) are labeled by the eigen- 
value of the quadratic Casimir, whose action on one-particle states is given by the de Sitter 
Laplace-Beltrami operator. Thus, one-particle states satisfying (2.3) have Casimir eigen- 
values —a{a + d). There is a redundancy in this description, as (2.3) is unchanged under 
the interchange 

a^-{a + d). (2.9) 
Unitary representations must have real Casimirs which requires 

(T = CT, or a + a = —d,. (2-10) 

This gives three distinct series of representations [38, 39] 

1. Principal series: representations with a = + ip, /) G M. 

2. Complementary series: representations with —d < cj < 0, o" G M. 

3. Discrete series: representations with cj > — |, cr G Z. 

If we restrict attention to the non-tachyonic cases ^^(M^ -|- ^K) = —a{a -|- d) > and use 
(2.9) we can restrict the range of a to 

-^<cT<0, 0- = -hi/O, p G M . (2.11) 

In particular, the "massless" scalar field with £^(M^ -|- ^i?) = corresponds to fj = 0; 
increasing £^(M^ -|- ^i?) moves a along the real line from to — (i/2; further increasing 
1'^{M'^ -\- ^R) moves a off the real line along a = — | -|- ip. 
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Finally, let us return to the equation of motion (2.3) and construct an explicit set of 
mode functions ip ^ for each a in the above range. To do so, let ('^^V^ be the Laplacian on 
S and note that the rescaled wave operator 



'a 



-d^ - (i(tanh t) dt + cosh^ t ('^^ 



(2.12) 



is the analytic continuation to dSo of the Laplacian (^)V^ on . As a result, solutions 
to (2.3) are analytically continued scalar spherical harmonics. Such spherical harmonics 
are reviewed in Appendix A.l and lead to a basis of solutions on dSu of the form 



(2.13) 



Here ^'^^Y^iyid) are scalar spherical harmonics on the spatial S'^ slices, while the time- 
dependent factor in (2.13) is given by [42] 



y^jit) = iVd,j(cosht)-(^-i)/2p- 



(i+{d-l)/2) 
f(d-l)/2 



(isinht) . 



(2.14) 



In (2.14), Py ^{x) is the Legendre function [43] and N(i„j is a normalization constant 



N, 



daj 



r(j - (j)r{j + a + d) 



1/2 



(2.15) 



Note that N^^^j is real for all values of cr (2.11). 
2.2 Boost matrix elements 

We are now ready to compute matrix elements of the finite boost B{X). In order to compare 
with the vector and tensor cases in section 3, it is useful to begin with the infinitesimal 
generator B. To this end, consider the de Sitter Killing field 



= (cos X; — tanh t sin x, 0) 
The corresponding boost generator 13 acts on ip^ as 



B,(po 



i£BP>a = i {cosxdt — tanh t sin x^x) Po 



(2.16) 



(2.17) 



Using the well-known properties of Legendre functions [43] one finds 



i£B^ ^ = h'^^'P ^,1+6 -^<p -f 1 



(2.18) 



where j + 1 and j — 1 are j with j replaced by j it 1 and all other unchanged. The 
coefficients b^^ are 



U + k + d- -k + l)(j - a){j + g + d) 
(2j + d-l){2j + d+l) 

^j + k + d- 2){j -k){j-l-a){j-l + a + d) 
(2j + d-l)(2j + d-3) 



1/2 



(2.19) 
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Note that 6 - = when j = k so that the coefficient of w ^ -, vanishes when u? , is not 
defined. As expected, the boost generator 

^ = -E«K/«M + ^J«J+i)' (2-20) 
1 

is Hermitian. 

Let B{X) = e^^^ be the operator which translates the field by an amount A along 
(2.16). -B(A) acts on scalar fields via the coordinate change a; — > x' given by 

sinh t' = cosh A sinh t + sinh A cosh t cos x > 
cosh t' cos X = sinh A sinh t + cosh A cosh t cos x > 

cosh i' sin x' = cosh t sin % . (2.21) 



We wish to compute the matrix elements 



BjJX) := (i?(A)^^.,^, 



'KG 



(2.22) 



of -B(A) between one-particle states associated with the mode functions (2.13). Due to 
rotational symmetry about the boost axis, the matrix elements (A) are diagonal in all 
angular momenta except the total angular momenta j, m; i.e. Bj^{X) oc Sj-^mi ■ ■ • ^jd-i,md-i ■ 
Furthermore, the matrix elements are symmetric in j <^ m and depend only on j, m, and 
the (mutual) second quantum number k = jd-i = w-d-i- See [38] for a group theoretic 
explanation of these facts. As a result, we may adopt a simpler notation for the matrix 
elements, denoting them by Bjmk{^) = Bmjk{^)- 

Using the explicit form (2.13) of the modes and the coordinate change (2.21), one 
can compute the inner product (2.22) directly. This lengthy calculation is performed in 
appendix B and yields 



C [{m-k)/2] [{j-k)/2] 

i?,mfe(A) =^J^,(coshA)^-'=(isinhA)/ Yl E Bf^ki^oshX)' 

{ a=0 b=0 



2b 



X 3F2 



2 



j+a+d+f 
2 



j+m+f+l 
2 



k — a — b: 



■ sinh A 



1 I f j+m+d+l+f . 

2 J ' 2 ' 



•. (2.23) 



Here 3-^2(05 b, c; d, e; z) is the generalized hypergeometric function ([44, 45]) the coefficients 
A'^^f^ and ^^^fc are defined in (B.20) and (B.21), and / is defined below (B.16). We 
conclude this subsection some comments on this result: 



• Due to some asymmetries in our treatment of the j,rh mode functions, the result 
(2.23) is not manifestly symmetric in j and m. However, for any given example one 
can readily verify that (2.23) is symmetric by, e.g., using Mathematica to expand 
(2.23) as a power series in sinh A. 
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It is well known that the boost matrix elements satisfy a number of recursion relations 
[25, 46]. While we will not use any such relations here, we note one of the most useful 
relations here in order to connect with the work of previous authors. 



(2.24) 



The result (2.23) is rather cumbersome in general, but simplifies greatly for two cases. 
The first is the case where the scalar field is conformally coupled, a = —{d — l)/2, 
for which (2.23) reduces to 

j+m+d—l 



sinh — I I cosh — 
2 \ 2 



X 2F1 



j + k + d— — k + 1; j — m + 1; — sinh' 



A 



(2.25) 



where in general we define 



Pjmk{o') 



■j-ni 



2J-™r(l + j - m) 



j + k + d- l,j-k + l,j-a,j + a + d,m + + 



d+l 



m + k + d— l,m — k + l,m — a,m + a + d,j + , j + 



d+l 



1/2 



(2.26) 



Here r[. ..] denotes a product (and quotient) of gamma functions; e.g T 



ah 
c 



^^T[<^^ • I^oth (2.25) and (2.26) assume j > m. A second case in which the matrix 
elements simplify is when m = k. In this case (2.23) becomes 



BjkkW = BkjkW 

= /3,fcfc(a)(sinhA)^-SFi 



j — cr j + a + d d+l 

- ; J + 



d + l 



sinh A 

(i+{<i-l)/2) 



/3,..(a)(sinhA)-'^-('^-^)/^P;^{+it-^^^ 

(2.27) 

Matrix elements for nearby Bj^k+n,k{'^) can then be found by successive use of the 
recursion relation (2.24). 

2.3 Convergence of group averaging 

We wish to use group averaging over the de Sitter group G = SOo{D,l) to construct 
physical states which satisfy Ma_b\^) = 0, where Mab are the generators of the group. 
It is useful to first find simple expressions for unitary transformations U{g) and the Haar 
measure dg of the de Sitter group. Under the Cartan decomposition of SOo{D, 1) we can 
decompose any group element into two SO{D) rotations and a boost [38]: 



U{g) = U{a)B{\)U{^) , 



(2.28) 
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where a and 7 are group elements of SO{D) C SO{D, 1) and B{X) is the unitary action in 
the given representation of a boost of rapidity A along some given Cartesian axis in M^j+i. 
We take the coordinate x to have values x = 0; on this axis. In a similar fashion, the 
Haar measure can be decomposed as [53] 

dg = da dj d\ (sinh A)'^ , (2.29) 

where again d = D — 1 and da,d'j are both the Haar measure on SO{D). The group 
averaging inner product (1.2) is then 



(^i|^2) = / dad^ dX {sinh Xfi^PijU {a) B{X)Uij)li^2) 

dX{smhXf{%bilVoB{X)Vol^/j2) , (2.30) 



where in the second line we used the fact that the projector Vq onto SO(D)-invariant states 
satisfies 



Vo = j da U{a). (2.31) 

Due to the projectors Vq, it is only necessary to compute (2.30) for rotationally invari- 
ant auxiliary states IV'1,2)) in which case the action of the projectors is trivial. In practice, 
it will often be useful to simply delete the projectors and to consider the integral 

POO 

1= dA(sinhA)'^(^ilS(A)lV'2) . (2.32) 

JO 

Because SO{D) is a compact group, the operators Vq do not affect the convergence prop- 
erties of the group averaging inner product. Examining the convergence properties of X is 
thus equivalent to examining those of (^'i|^'2)- 

For A^-particle states in which each particle occupies a definite mode, the integral I 
may be written 

/•oo 

1= / dA(sinhA)'^(ii,...,jjvl5(A)lmi,...,m^) 
Jo 

POO 

= dX{smhXrBj^^.^{X)---B^j^^.JX), (2.33) 



where the product of matrix elements B^ ^ • • • B^ ^ in the last line has been sym- 
metrized over the indices mi • • • rh]\f. For large A the boost matrix elements have the form 



BjmkW = Co,mfc(e^)"-' 1 + Cij^kie-''^) + 0{e 



-2A^ I r)f^-4:X\ 



+2?0imfc(e^)~''^"^'' [1 + Vijmde~'l + 0(e-^^)J (2.34) 

where Ci are coefficients that do not depend on A and Vi = Ci{a ^ —{a + d)). Therefore, 
the least convergent part of the integrand in (2.33) has the asymptotic behavior ~ exp[A(ci-|- 
Na — ki)]; the integral converges for 

Re^d + Na -Y^k^ < . (2.35) 
This has the following results for scalar representations: 
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• Principal series: The integral (2.33) converges absolutely for > A'^conv = 2. 

• Complementary and Discrete series with a ^ 0: The integral (2.33) converges 
absolutely for N > Nconv = —d/a. This can in principle be a very large number for 
a near 0. 

• Massless scalars: The matrix elements Bjmk{^) above do indeed form a represen- 
tation of the de Sitter group for a = 0. However, this representation does not include 
the full physics of the zero mode as this representation describes only states with 
vanishing zero-mode momentum; see [20]. 

In all cases the integral may converge for N < Nconv for sufficiently large momenta ki. 
2.4 Convergence for more general states 

The convergence criteria listed above have been derived for states in which each particle is 
described by a single mode in our basis. One would also like to examine more general states 
which allow the particles to be in a superposition of modes. Finite superpositions clearly 
do not alter the above criteria. However, understanding the convergence of infinite super- 
positions requires some control over the group averaging norms of the states . . . , jat) 
as a function of the angular momenta. To see this, let us consider an arbitrary smooth 
A^-particle state fV'Tv) which may written as an infinite sum 

1^^)= E ^Ji,...Jivl^i'---'^") (2.36) 

with coefficients cj^ which decay with all ji faster than any polynomial. The group 
averaging inner product of two such smooth states will be finite so long as the inner product 
of the basis states grows no faster than a polynomial in momenta ji. 

In principle one can compute the group averaging inner product (1.2) exactly using the 
finite boost matrix elements (2.23) and SO{D) Clebsch-Gordan coefficients. Unfortunately, 
it is rather difficult to extract the angular momentum dependence of this result; indeed, 
even if one omits the SO{D) Clebsch-Gordan coefficients as in the previous section, it is 
difficult to extract the angular momentum dependence of the finite boost matrix elements 
BjmkW- In order to proceed, in this section we derive a simple bound on the finite boost 
matrix elements and use this bound to investigate the dependence of the group averaging 
inner product on angular momenta j, m, A;. The bound's benefit is that it's dependence 
on angular momenta is clear; it's weakness is that it has a weaker fall-off with respect to 
rapidity A than the actual matrix elements. As a result, larger numbers of particles will be 
needed to guarantee convergence for general smooth wavefunctions. 

The calculation of the bound on Bjmk{^) will be somewhat technical. We therefore 
begin by summarizing our results, and then follow with the derivation. We are able to 
show that the matrix elements \Bjmk{^)\ ^ /(i) k)g{X), where /(j, m, k) contains all the 
angular momenta dependence and g{X) contains all the rapidity dependence. The function 
f{j, m, k) grows at large momenta no faster than a polynomial of order 0{{j'm)^'^^^^^'^k~'^). 
The function g{\) has asymptotic behavior at large A ^ 1 given by 
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j e ^ for (/c- Re cj) > 1 
9(A)~ for(/fc-Rea)<l ^^"^^^ 

Recall that the exact expression for the matrix elements decays as ~ exp[X{a — k)]. Bound- 
ing Bjmki^) in this way thus looses the additional suppression that occurs for large mass 
and large angular momenta k. However, it allows us to conclude the following for arbitrary 
smooth states: 



• Principal series: For d > 2 principal series representations always satisfy Recr < 
— 1. Thus for this series our bound on the matrix elements decays with A like ~ e~^. 
It follows that for N > d group averaging converges for any smooth states. 

• Complementary series with o" < — 1 or A; > 0: For d > 2 the complementary 
series contains representations with a < —1. The conclusion is as for the principal 
series: for N > d group averaging converges for any smooth states. 

• Complementary series with a > —1 and k = 0: For small a, i.e. light mass, and 
A; = our bound decays with A like ~ e'^'^, so that group averaging converges for any 
smooth states with > —d/a particles. This agrees with the threshold Nconv found 
in the previous section for finite-superposition states. 



We will discuss the status of states with smaller numbers of particles in section 4. 

In the remainder of this section we derive the bound on the finite boost matrix elements 
BjmkW discussed above. From the definition of the boost matrix elements (2.22) we have 



'KG 



-2i 



(2.38) 



t=0 



In the second line we have integrated over the S'^~^ leaving only the t- and x-dependant 
parts of the modes, and have also used the fact that when evaluated on positive frequency 
states the two terms in the Klein-Gordan inner product are equal. The function ^'^^yjkix) 
is defined in (A. 2); the function yaj{t) is defined in (2.14). To simplify the remaining 
integral over x we will bound both '^'^^yjkix) and j><xj(i) 

in turn. 

We begin by writing ^"^^yjkix) iii terms of Gegenbauer polynomials C^(x): 



id) 



yMx) 



T{2k + d-l) 

2fc+(d-l)/2r(A;+ I) 



(2j + d- i)r(j -k + i) 

T{j + k + d- 1) 



1/2 



ismx)'c':^^'- 



l)/2, 



cosx)- 
(2.39) 
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We may bound the size of ^'^^yjkix) via the addition theorem of Gegenbauer polynomials 
[38], which may be written 



Cj(cos 9 cos X + sin 9 sin x cos ip) 
j 

= ^22^(2p + 2A;- l)r 

k=0 



p + m,p + m,2p — — m 

p, p,2p + j + k 



X (sinxsin0)'=Cf+ncosx)C^:!"fc(cos0)C^^'/'(cosV'). (2.40) 



If we set -0 = and 9 = x then each term is positive and we have that 



(sinx)'^C'=+,^(cosx) 



< ' r 

- 22fc(2p + 2A;- 1) 



p, p,2p + j + k 

p + m,p + m,2p — 1,1 + j — m 



(2.41) 

valid for < k < j, p > 1/2. From this expression we can derive two useful bounds on 
('^)3^jfc(x). First, setting k ^ and p ^ {d - l)/2 in (2.41) we find 



|('')3^,,.=o(x)l < 



(2j + d- l)r(i + d-l) 
2d-^T{l + j)T (If 



1/2 



(2.42) 



This expression has technically been derived only for d > 2 (recall we must have p > 1/2 
in (2.41)); however, (2.42) turns out to be valid for d = 2 as well. Second, for the case of 
A; 7^ it is useful to have a bound that includes some remnant x-dependence. For this case 
we let p ^ (d - l)/2 + 1, A; ^ A; - 1, j ^ j - 1 in (2.41) and find 



(2j + d-l) . 
^{2k + d - 2)d 



k,j + d,^ 
j,k + d-l, 



1/2 



(sinx). (2.43) 



This expression is valid for d > 2 as well. 

Next we examine the time-dependant part of the modes 3^o-j(i)- First, from (2.14) we 
compute 



dtyamit) 



(2.44) 



To bound yaj{t') we must bound the Legendre function P^^^^'^^'^^y^^'^\ismh.t). We may 
define the Legendre function through the integral representation 



l)A'/2 



da;(z+coshx)''"''~^(sinhx)^''+^ Re(-/i) > Re{u) > -1 



2u+l 



for Re(— //) > Re(z^) > — 1; using this we bound 



\y^j{t)\ < 



-^^"^1 - |2.+(d-i)/2r [^+d±i,j_^j 



(2.45) 
(2.46) 
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where 



-(i+Re cr+d)/2 



(sinh X 



.2Rea+d 



(2.47) 



For a > —d/2 (corresponding to the complementary series) it is easy to bound / by 
(sinhx)^'^'*"'^ < (cosh x) (sinh x)^'^"'"'^"^ and using the identity: 



poo 

/ (cosh x) (sinh x)^"~^ (1 + zsinh^ x) 
Jo 



2 ^\-a-b 



a, b 
a + b 



for Rea > 0, Re6 > 0, z > 0, yielding 



\yaj{t)\ 



N. 



daj 



2a+{d+l)/2 



(cosht)'^, cr > 



(2.48) 



(2.49) 



For the case Retr = —d/2 (corresponding to the principal series) we must bound I differ- 
ently; here 



/ = (cosht) 



-d/2 



Now note that 



so that 



dx (l + sech^tsinh^ x) 
1 1 



-(i+i)/2 



(l + sech isinh x) < — \- -e 
^ ^24 



2{x-t) 



(2.50) 
(2.51) 



/ < (cosh t) 



-d/2 



dx 







1 1 



2{x-t) 



-1/2 



2(cosht)-^-'^/2arcsh(\/2e*), (2.52) 



and we have the bound 



\yaj{t)\ < 



:Y— (cosht)-'^/2arcsh(\/2e*), (Rea = -d/2) . 



(2.53) 



We conclude by combining (2.49) and (2.53); because we are primarily interested in the 
behavior of Bjmk{^) at large j » 1, we take the liberty of assuming (j — Recr) > 3 which 
allows us to simplify the expression slightly: 



\yajit)\ < 



1 



|r(a + ^; 



j"0-+l j+a+d+1 

2 ' 2 
j-T j+er+d 
2 ' 2 



1/2 



(cosh t) 



Re <T 



ir {a + for a > -d/2 
/2 arcsh( e*) for a = -d/2 + ip ' 



(2.54) 



Note that in this expression the t- and j- dependence separate. 

We now insert our bounds (2.42), (2.43), and (2.54) into our expression for the finite 
boost matrix elements (2.38). From the coordinate transformation (2.21) we see that at 
t = we have (cosht')^ = 1 + sinh^ Acos^ x and that t' < A; using these facts one may 
evaluate the integral (2.38), finding the bounds 



\Bj,m,k=oW\ < fU,m,k = 0)go{X), |-Bj,m,fc>o(A)l < fij^rn^k > 0)5ri(A), 



(2.55) 
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where 



f{j,m,k = 0) : = 



|r(a + ^) 



j — cr+l j+a+d+l m—a+1 m+a+d+1 

2 ' 2 ' 2 ' 2 
j — cr j+a+d m—a m+a+d 

2 ' 2 ' 2 ' 2 



1/2 



j + d — l,m + d — 1 
1 + j, 1 + m 



1/2 



f{j,m,k> 0) : = 



|r(a + ^) 



i-o-+l j+fJ+rf+l m-cr+l m+o-+d+l 

2 ' 2 ' 2 ' 2 
j — c j+a+d m—a m+a+d 
2 ' 2 ' 2 ' 2 



r(A;)A/(j + 



2 / 



[m + 



d-l^ 



(d + l){2k + d- 2)T{k + d-l) 



j + d,m + d 
j,m 



, (2.56) 

1/2 



1/2 



, (2.57) 



go{X) := (cosh A) 



5i(A) 



Rc a 



2F1 



Re (T d d + 1 



; tanh^ A 



((J + ^) for a > -d/2 



I \/2 arcsh(\/2 e^) for a = -d/2 + ip 
(cosh A) ^^'^-^ 



2^1 



1-Recj d + 2 d + 3 ,2, 
, — - — ; — - — ; tanh A 



(a + 



for a > -d/2 



^/2 arcsh(^/2 e^) for a = -d/2 + ip' 



(2.58) 



(2.59) 



Applying Stirhng's approximation to the gamma functions confirms that /(j, m, k) is 
bounded by a polynomial in j as desired. The asymptotics of g{X) state in (2.37) fol- 
low from the well-known asymptotics of the hypergeometric function [43, 44, 45]. 



3. Higher Spin 

The techniques of section 2 can also be used to study the convergence of group averaging 
for linear vector and tensor gauge fields. In fact, with the appropriate choice of gauge and 
in the appropriate basis, the matrix elements of the infinitesimal boost generator B will 
turn out to have precisely the same form as in the scalar case. As a result, we can use our 
results from section 2 to evaluate the matrix elements of -B(A) between vector and tensor 
one-particle states and to quickly arrive at the desired result. 

3.1 Massless Vector Fields 

Massless vector fields obey the Maxwell equation of motion 

OA^ - ^A^ - V^V'A^ = , (3.1) 
which is of course invariant under the gauge transformation 

A^^A^ + iVf,^ , (3.2) 
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for any scalar ip. Note that V^^^ — > + so choosing Lp = I □ V^^^ brings 

us to a gauge in which Au satisfies 



□ 



0, V^^„, = . 



(3.3) 



There remains some residual gauge symmetry which we will discuss below. As in the scalar 
case it is convenient to define a quantity ai related to the eigenvalue of the wave operator. 
For vector fields the natural choice [42] is 

-cri(cri + d) + 1" 



UA,. 



(3.4) 



so that (Ti = — 1 for the massless case. 

Transverse vector fields A^ satisfying (3.4) in dSo are the analytic continuation of 
transverse vector spherical harmonics on (reviewed in Appendix A. 2). There are D — 1 
sets of independent modes, corresponding to the fact that each set of modes has non- 
vanishing components only for indices tangent to the sub-sphere 5^, 5^, ... , S^. We label 



the sets of modes A 



(2;i) 



^if where J is as in (2.5). 



As mentioned above, there is some residual gauge freedom in (3.3); namely, there 
remains the freedom A^^ — > -|- £V^(/9o for ipo satisfying the Klein-Gordon equation with 
fj = 0. In appendix A. 3 we use this freedom to eliminate the modes A^ 



{D;j) 



Since Af 



vanishes for all other modes, we have fixed temporal gauge At 
are 



0. The remaining modes 



AiP-J) = ^(3-rf)/2(^oght) %,j{t) ('^)yi^'^')(0,). (3.5) 

Here a is a spatial index, (3 = 2, . . . ,d denotes the set of modes, and '-^^ yi^'"'^(^^d) are vector 
harmonics on S"^ (see Appendix A. 2). Note that such harmonics exist only for j > k > 1. 

As with our scalar modes, these modes can be shown to satisfy a de Sitter invariant 
positive frequency condition. In addition, they have been normalized with respect to the 
inner product 



KG 



(3.6) 



which agrees with the symplectic structure of the Maxwell field when At = as above. 
We therefore expand the field in creation and annihilation operators a^^^'^^ and a^^'^^ 

as 

M^) = E E U'^^^^H^) + a^^^^nfh-)] ■ (3.7) 

It follows from (3.6) that the creation and annihilation operators have the usual commu- 
tators 



«t(«;i) „t(/3;-) 



0, 



and allow us to define an auxiliary space vacuum state which satisfies 

a(^;j")|0)=0 V/3, J. 
General states in Waux are linear combinations of the states a'^^'-'^^ • • • a(^"'-^")|0). 



(3.8) 
(3.9) 
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3.2 Vector Boost Matrix Elements 

We now compute the matrix elements of the boost generator B that generates boosts along 
B'^ (2.16). The action of B is given by 



B,A, 



i£BA^-P.G., (3.10) 



where P.G. denotes the subtraction of a "pure gauge" term required to maintain the above 
gauge conditions. In particular, this subtraction is required to maintain the temporal gauge 
condition At = which breaks manifest de Sitter invariance. It is convenient denote the 
action of a boost on the modes by 6b so that 

B, ^m] = E E «^^''^'5b4^''^ + h.c. (3.11) 

(3 ■ 

In general, this action takes the form 

6bAI^'^^ = i£BAf^'^ - clV^if^j , (3.12) 

where c is a constant. 

Let us first consider the modes ^4^^'''^ for (3 ^ d. Since these modes vanish in the 
subspace of the co-tangent space associated with (i.e., a[^'''^ = and A^^'''^ =0), we 
also have 

£BAf'^^ = £bA^^'^^ = . (3.13) 

Thus, the Lie derivative preserves temporal gauge when acting on such modes and there 
is no pure gauge contribution; the action 5b is given simply by i£B- After a short 
calculation one finds 

6bA(^'^^ = b+ .yl(f'J+i) + b- .^|f'^"-i) , (3.14) 

where the coefficients are identical to those for scalars (2.19) with a = —1. These modes 
"transform like scalars" in the sense that they have boost matrix elements identical to 
scalars with a = —1 (as one might expect for massless vector fields). 

It remains only to consider the modes Once again, the details are presented in 

appendix A. 3 where it is shown that SbA^^'''^ takes the form (3.12) with 



(cosh tdt + (d-l) sinh t) ^^Ai^'^^ 

cosh t ^ 



(3.15) 



It is now straightforward to compute matrix elements of B for the A^f^''^^ modes. A key 
point is that the Lie algebra requires SbA^i^'"^^ to be a linear combination of j4|f ""^^ and 
A^f''' ^\ so that we need only determine the coefficients. This is most easily done by 
considering the single component Ai^'''\ After some tedious algebra one finds 

5s4^'^') = b\^A[^'^+^^ + 6:^.4^'^'-!) , (3.16) 
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where the coefficients are identical to those for scalars (2.19) with a = —2. 

For each set of modes above, we found that B has precisely the same matrix elements 
as for scalar fields with either cr = — lorcj = — 2. As a result, the matrix elements of finite 
boosts B[X) can be read directly from ^ (2.23) and the group-averaging properties follow 
immediately. The results for vector fields are thus: 

• A^^'^'^^: These modes act like scalars with a = —2. Thus the group averaging expres- 
sion I (2.32) for A^-particle states built from a finite number of such wavefunctions 
converges when d-2N - J2i h < 0. Since the minimum value of k for these modes 
is A; = 1, this generally requires > Nconv = d/3. 

• A^^^'-') for /3 7^ d: These modes act like scalars with a = —1. Thus, for A-particle 
states built from a finite number of such wavefunctions, I converges when d — N — 

ki < 0. Since the minimum value of A: is A; = 1, convergence generally requires 
N > Nconv = d/2. 

• Both: For A-particle states built from a finite number of wavefunctions of either 
type, I again converges when N > Nconv = d/2. It is also clear that I may converge 
for A^ < Nconv if the particles have high enough angular momentum k and/or mostly 
have (3 = d. 

• Arbitrary smooth states: An arbitrary smooth wavefunction may be written as 
an infinite superposition of the above modes, each of which has a < —1. From 
our discussion in section 2.4 we conclude that group averaging converges for general 
smooth states when N > d. 

3.3 Gravitons 

We now address group averaging for free gravitons, following the same recipe as in the 
preceding sections^. After quantizing metric perturbations to create the auxiliary Hilbert 
space "Haux, we compute the infinitesimal boost generators for gravitons and find that they 
again agree with expressions found in the scalar case. Thus, as in the case of vector fields, 
the results of section 2 allow us to quickly arrive at conditions under which group averaging 
will converge. 

We wish to consider linearized perturbations /i^j^ of the background de Sitter metric 
g^i^. Starting from the Einstein-Hilbert action S = J d^Xy/—g{R — 2A), one finds the 
linearized equation of motion 

= DV + V,,V^/i - V^V^hxu - V^V^hx^, - ^ V 

+ (^V"V^/i«^ -ah+ ^-j^h^ g^, , (3.17) 

^Note that for D = 4 (i.e., d — the matrix elements with a = —1 and a = —2 are identical as a result 
of the invariance under cr — (a + d), and are identical to the matrix elements of a conformally coupled 
scalar field. They therefore reduce to the simplified expression (2.25). 

^Group averaging of gravitons for d — 3 was considered in [9]. For this special case, one can obtain more 
detailed results than we derive below. 
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where h = g^'^h^u and once again □ = V^V^ is the wave operator on de Sitter. We impose 
transverse traceless gauge 

Vxh^^ = 0, /i = 0, (3.18) 
so that the equations of motion become 

° - f ) V = , V^/iam = /i = . (3.19) 

Not surprisingly, mode solutions to (3.19) are transverse traceless tensor spherical 
harmonics on analytically continued to dSo- Transverse traceless tensor harmonics are 
reviewed in Appendix A. 4. After analytic continuation, the general such harmonics satisfy 

"''^t"'^' ] V . (3.20) 

in addition to the transverse-traceless condition. Solutions to our equation of motion (3.19) 
are those with o"2 = 0. There are {D + 1)(-D — 2)/2 sets of independent modes with non- 
vanishing components tangent to the sub-spheres 5^, . . . , . We label the sets of modes 
where a, /3 = 2, 3, . . . , d, D. Individual modes are then labeled by momenta j. 
Recall that for free vector fields imposing transverse gauge still allowed the residual 
gauge symmetry — > + clV ^ipo^ and that used this symmetry to set At = (temporal 
gauge). A similar story holds for gravitons. The equations of motion (3.19) are invariant 
under the gauge transformation h^j_^ h^y + ^V(^Aj^) so long as obeys 

° + 4 = ' ^''^/^ = • (3-21) 



Following the same steps as in the vector case (see Appendix A. 3) one can show that the 
tensor modes h^J^"\ a = 2, . . . , D are in fact proportional to V(^A^), and can therefore be 
set to zero, fixing the gauge completely. Since these are the only modes with non-vanishing 
t component, this amounts to fixing htfj, = 0. 

The D(D — 3)/2 remaining sets of modes are 

'hfi — ^ ' 

h^T^^ = £^'-'y\cosht)'%j{t)^'^^Tif'^\n,), a, P = 2,...,d, (3.22) 

where Wrf,"^'^')(17rf) are tensor harmonics on S (see Appendix A. 4). Note that rank two 
tensor harmonics exist only for jii > jd-i > ••• > J2 ^ \ji\ ^ 2. As in the scalar and 
vector cases considered before, these modes satisfy a de Sitter invariant positive frequency 
condition. The modes have been normalized with respect to the inner product 

(/,{!), := I d^^h^^}vlWr . (3.23) 
We can therefore expand h^^ in creation and annihilation operators 



" /3 7 



(3.24) 



and complete the canonical quantization of /i^i/ much as was done above for 99 and A^] see 
[47, 48, 49, 50, 51] for useful references. 
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3.4 Graviton boost matrix elements 



Finally, we compute the finite boost matrix elements for gravitons. The calculation is 
essentially the same as for free vector fields, so we will simply summarize the results. 
The boost generator B acts on the graviton field as 

= i£Bh^u -P.G. , (3.25) 

where once again P.G. denotes a pure gauge term that must be subtracted in order to 
preserve temporal gauge ht^ = 0. The pure gauge term must be proportional to ^V(^yli,) 
where satisfies (3.21). The calculation of this pure gauge term is analogous to the 
calculation done for free vector fields presented in appendix A. 3. One finds that all modes 
transform like 




with the following values of a: 

• hfj^i^^^'^'^^ : these modes transform with a value of cr = —2. Thus, for A^-particle 
states built from a finite number of such wavefunctions , I (2.32) converges when 
d — 2N — ki < 0. Since the minimum value of k for tensor modes is A; = 2, this 
requires N > Nconv = d/A. These are the most convergent modes. 

• h^^^'^"'''\ a 7^ d: these modes transform with a value of o" = —1. Thus, for A^- 
particle states built from a finite number of such wavefunctions, I converges for 
N > Nconv = d/3. 

• h^i,^"^'-''^ a,/3 7^ d: these modes transform with a value of cj = 0. Thus, for N- 
particle states built from a finite number of such wavefunctions, 2 converges for 
N > Nconv = Nconv = d/2. These are the least convergent modes. 

• Arbitrary smooth states: For all modes the quantity a — k < —1, so following the 
discussion in section 2.4 we conclude that group averaging will converge for arbitrary 
smooth states when N > d. 

As in the vector case, modes orthogonal to are the least convergent while those 
with non-zero overlap are more convergent. For finite superposition states, the number of 
particles required for convergence is set by the least convergent modes to be A > Nconv = 
d/2. It is again useful to note that for D = 4 (i.e., d = 3) the finite boost matrix elements 
are the same for a = —2 and a = —1 and are given by the simplified expression (2.25). 

4. Discussion 

We have studied the general behavior of matrix elements in a standard basis for free 
massive scalars, abelian vector gauge fields, and linearized gravitons in any dimension d. 
We have shown that group-averaging converges for states with A-particle wavefunctions 
given by de Sitter harmonics - as well as finite superpositions of such basis states - when 



B,h 
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all components of the quantum state contain N > Nconv particles. For principal series 
scalars Nconv = 2 (so we require > 3), for scalars in the complementary series and in 
the discrete series with 7^ we have Nconv = d/\a\, while for both vector gauge fields 
and linearized gravitons Nconv = d/2. 

We also considered states with general smooth wavefunctions. Such states are in 
general infinite superpositions of the above basis states. Due to the difficulty in computing 
the group averaging inner product exactly, our results for this case relied on finding a bound. 
The key step here to bounding the finite boost matrix elements BjmkW as discussed in 
section 2.4. Our bound grows like a polynomial in total angular momenta but decays 
with rapidity at a weaker rate than the actual matrix elements. As a result, for principal- 
series scalars, linearized gauge fields, and linearized gravitons we can only demonstrate 
convergence for general smooth wavefunctions with N > d particles, a somewhat larger 
number than for finite superpositions of our de Sitter harmonics. We have by no means 
shown that such a high particle content is necessary. In fact, as discussed below, we suspect 
that the values of Nconv stated above for finite superpositions should also suffice for states 
with arbitrary smooth wavefunctions. 

Our bound on the finite boost matrix elements grows with angular momenta like a 
polynomial of order 0{{jm)^'^'^^^/'^), suggesting that the group averaging inner product 
should likewise grow with angular momenta. However, in cases where the inner product 
has been computed exactly (or numerically) it has been found that the inner product 
actually decreases with angular momenta. For example, Higuchi [9, 19] computed the 
group averaging inner product exactly for massless 1 + 1 scalars^, conformally-invariant 
scalar fields in 3 + 1 dimensions, and 3 + 1 massless tensor fields. All of these cases have 
a = —{d — l)/2, i.e. the same as for a conformally coupled scalar field. For this special 
value of a it is possible to define a set of states in Waux which map under group averaging 
to a complete orthonormal basis in Wphys for states with N > d/a = 2d/{d— 1) (and in the 
case of d = 1 Higuchi finds > 1); see also the discussion in [20]. Higuchi's results show 
that the group averaging norm decreases with total angular momentum, so that group 
averaging in fact converges for any normalizable state in Waux with > 2d/ {d — 1) (or 
iV > 1 for d = 1). 

The difference between Higuchi's results and our own might be attributed to i) the 
possibility that the solvable cases a = —{d — l)/2 may be exceptional, ii) some weakness 
of our bound on the matrix elements Bj^y^i^^X), or iii) the fact that we have not considered 
the effect of the SO{D) projection operators Vq- Numerical results suggest that (iii) is 
the most important effect. Consider for example a free 2 + 1 scalar field with a = —3/2 
(at the border of the principal and complementary series). The quantity I (the group 
averaging inner product without the 5*0(3) projection operators) from (2.32) is shown in 
figure 1 (left) for a particular sequence of normalized states. It is clear that 2 increases 
with j (though notably not as rapidly as our bound). In contrast, one can also consider 2 
for a family of rotationally-invariant 3-particle states of the form |2 singlet, j) (8> |1 singlet), 
where |2 singlet, j) is the unique normalized rotationally invariant 2-particle state in which 

■^See [20] for a treatment of the zero-mode. 
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Figure 1: Left: The group-averaging expression X = / (iA(sinh A)''(i?jjfe(A))^. Shown (from top 
to bottom) are the cases fc = 0, 2, 5, 10. Right: The group-averaging norms of rotationally-invariant 
states. The dashed hne is the best-fit curve ~ 3~^^'^ ■ For both figures iV = 3, d = 2, cr = —3/2 



each particle has total angular momentum j and |1 singlet) is the normalized rotationally- 
invariant one-particle state. Because such states are already rotationally invariant, X coin- 
cides with their group averaging norm. As shown in figure 1 (right), these norms decrease 
with J, and in particular appear consistent with a power law of the form ~ i^^^"^ ■ From 
this we conclude that including the 5*0(3) projection matrices is essential to obtaining the 
correct angular momentum dependence of the group averaging norm. If indeed the group 
averaging inner product of normalized rotationally-invariant states decreases with j, then 
the above-stated values of Nconv will guarantee that group averaging converges for general 
normalizeable states. 

We conclude by discussing the status of group averaging for states with N < Nconv 
Let us begin with the case = 0. In all cases addressed above, the field admits a de Sitter- 
invariant vacuum. It is clear that the group averaging norm diverges for this state, since 
(1.1) would integrate {0\U{g)\0) = 1 over the entire non-compact de Sitter group. On the 
other hand, since |0) is already de Sitter invariant, there is no need to group average it to 
define a physical state. This case is directly analogous to examples studied in [52] and the 
same conclusion holds: If O is de Sitter invariant, then so is OlO). But the only normalizable 
such state is the vacuum itself. Thus, under the action of de Sitter-invariant observables, 
one finds that the state |0) is superselected from all other states. Since group averaging 
is a technique for taking the action of this observable algebra on Ti^ux and constructing 
a representation on de Sitter-invariant states, it is clear that the vacuum must be treated 
separately. See [52, 14, 53, 22, 23, 24] for more discussion and further examples of this 
phenomenon. 

Another interesting point is that group averaging never converges for 1-particle states. 
One may see this from the results above by noting that all S'0(D)-invariant 1-particle states 
have j = and also noting that neither vector gauge particles nor linearized gravitons 
exist for D = 2. This lack of convergence for A^ = 1 is a natural consequence of the 
linearization-stability constraints, as they suggest that no physical 1-particle states should 
exist. In particular, at the classical level such constraints forbid the existence of 1-particle 
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spacetimes. By this we mean spacetimes which differ from de Sitter space only by a localized 
particle-like excitation carrying positive energy in the sense of the (timelike) weak energy 
condition. The constraints state that the de Sitter charges Q[£,] = f-^ n"'£,^Tab must vanish, 
while for a localized distribution of positive energy one may always choose to be timelike 
and future-pointing in the region where Tab is non-vanishing. There is nothing to 'balance' 
this charge on the opposite side of the Killing horizon (where becomes past-pointing 
or spacelike), so the constraints cannot be satisfied. 

Though group-averaging of our 2-particle basis sates converges for the physically in- 
teresting cases of D = 4 gravitons and vector fields, this is not generally true in higher 
dimensions or for massive scalars. Now, classical 2-particle spacetimes do exist, but they 
are highly constrained as the two particles must 'balance' exactly. In fact, in the test- 
particle limit the worldline of either particle must be mapped to the other under the de 
Sitter antipodal map. As a result, these worldlines preserve an Mx SO{D — 1) subgroup of 
the de Sitter group SO{D, 1). The M factor is associated with a non-compact boost gener- 
ator. Thus, at least in the classical limit, the behavior of two-particle states is much like 
that of the vacuum in that they become invariant under a non-compact symmetry group. 
This suggests that 2-particle states may also form a superselection sector under the observ- 
able algebra, and that some renormalized version of group averaging (see [14, 53]) might 
converge for such states. While dispersion of the wavefunction breaks the above symmetry 
in a quantum 2-particle state, it may be that this effect is generally small enough for the 
above analogy to hold. 

Cases where group averaging fails for A^-particle states with N > 2 may be similar. 
We note that this occurs only for scalars in the complementary and discrete series and for 
our massless higher spin fields. None of these cases have a classical 1-particle limit since 
their Compton wavelength is the de Sitter scale or larger. 

Though group averaging does not generally converge for N < Nconv, one may use the 
fe-dependence of (2.35) to show that group averaging does converge for certain carefully 
constructed such states. The trick is to build a state for which no particle has k = 0. 
We do not attempt to classify such states, but merely present a simple example. Recall 
that for all d the tensor product of two spin j representations of SO{d + 1) contains an 
SO{d+ l)-invariant singlet. In our basis, most components of this singlet state have k ^ 
for both particles, though one component has k = for both particles. However, we can 
use this 2-particle singlet to build several distinct 4-particle singlets. Let |(12)(34)) be 
the tensor product of the 2-particle singlet state formed from particles 1 and 2 with the 
2-particle singlet state formed from particles 3 and 4, and similarly for |(13)(24)). Then 
|(12)(34)) — |(13)(24)) is a rotationally-invariant state in which every component has k ^ 
for at least two of the four particles; the component with k = for all particles cancels 
between the two terms. In this sense we can always find SO(D)-invariant states in which 
X^j ki > N4/2, where is the largest integer less than or equal to N which is divisible by 
4. We can therefore find states where group-averaging converges whenever > 2d, even 
when a is very close to zero. 
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A. Symmetric tensor spherical harmonics of rank r < 2 



In this section we briefly review symmetric tensor spherical harmonics (STSHs) on the 
d-sphere S"^. Tensor harmonics of general rank are discussed in [42]. We also provide a 
more detailed discussion of spherical harmonics of rank r < 2. 

We use coordinates = {9d-i, . ..,9i,4>) on S"^. The metric on is simply dQ,f = 
d(j)'^, and on S'^ is 

dnl = dOl^^ + sin^ Od-idnl^^ 

= dej_^ + sin^ ed^i{dej^2 + sln^ Od^2{- • • + sin^ 92{dej + sin^ Oi d<j?) ...)). (A.l) 

A.l Scalar STSHs 

Scalar spherical harmonics on 5*^ are denoted ^'^^Y^ip.d)- They are labeled by d integers 
j = jd, • • • , ji with jd > jd-i > • • • > i2 > liil- The harmonics satisfy 

('^)v2('^)y.(0rf) = -jdUd + d- i)(^)yj(o<i) , (A.2) 

where ('^^V^ is the scalar Laplace-Bertrami operator on S^. For d = 1 this is^-'^^V^ = 9^ 
and for d > 1 we have 

^"^V^ = dl_^ + {d-l){coted-i)de,_, + ^^^ffj'y , d>l. (A.3) 



The solutions to (A.2) are then 

-iji<l> 



where 



/27r 



.n=2 



d>i, (A.4) 



(2j + n- l)r(j +k + n-l) 



1/2 



(sin ( 



2T{j -k + 1) 

Here Py^{x) is the associated Legendre function [43]. The function ^'^^yjkiO) is normalized 
such that 

(A.6) 



de{sme)''-'^''^yjk{e)(''^yfk{9) 



for d > 1. As a result, the scalar harmonics are normalized with respect to the inner 
product 

{Yj, Yj,)l2 := J dQd Y.{Vld)Y^Xnd) = Sjj, . (A.7) 
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A.2 Vector STSHs 



Following [42] we create vector STSHs on S"^ out of scalar and vector STSHs on S'^~^. It 
is convenient to denote the metric, covariant derivative, and Laplace operator on a sphere 
gn Mj^j^^ ^^"^Da, and (")V^. We also introduce the notation x = ^d-i- Using our 
coordinates (A.l) the line element on S"^ is 

ds'^ = ^'^'^-fabdx'^dx^ = dx^ + (sin xf ^'^'^hijdx'dx^ . (A.8) 

We let indices a, b range over S'^ while indices i,j will have range only over the 5"^~^. 
Vector STSHs satisfy 

id)j:,aid)y^^Q^ ('^)v2('^)K = [-j-(j + (i-l) + l]('^Va (A.9) 

and depend on d angular momenta j that satisfy jd > jd-i > • • • > j2 > \ji\- As in the 
main text we set j = jd and k = jd-i- There are d — 1 sets of orthogonal solutions to 
these equations, one set each with non-vanishing components tangent to the sub-sphere 
52, . . . , 5^. The sets are labeled ('^)y(2), My(3)^ . . . ^ {d)y{d) ^ ^nd individual harmonics 
are also labeled by their momenta j. These harmonics are given by 



^''W^^'^\nd) 



k{k + d-2) 
{j + l){j + d-2) 

1 



1/2 



(sinx)-' ^'"^y^kix) ^'"-'^Y^.i^d-i) , (A.IO) 
1/2 



{j + + d-2)k{k + d-2) 



X sm X 




(^x + {d-2) cot x) ^'^y^x) ^'-'^ A ^"-'^YMd-i) ,(A.ll) 



(nd) 



a ^ d 

sinx)^'^yMx)^'-'^vt''\nd^ 



a ^ d 



(A.12) 
(A.13) 



Since (A.IO) vanishes for A; = 0, we see that non-trivial harmonics exist only for j > k > 1. 
The vector harmonics are normalized with respect to the inner product 



(A.14) 



A. 3 Residual Gauge Transformations 



This appendix deals with certain issues related to the residual gauge transformations noted 
in section 3.1 preserving the gauge condition V^^^ = 0. Such transformations take the 
form — > Afj^ + i'Vfj^ipo, where (po satisfies the Klein-Gordon equation with a = 0. These 

transformations can be used to entirely eliminate the modes and also have implica- 



tions for the modes A 



We address each in turn. 



The modes A 



(D-j). 



The above residual gauge transformations can be used to eliminate the modes j4 
The point is simply that ^^^'■'^ oc V^c^g j, so that these modes are pure gauge. 
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We begin by introducing some new tools. Note that the following Casimir-raising and 
lowering operators 



act on yaj{t) as 



Lcr := i{coshtdt — crsinht), 

Ra := i{cosh.tdt + (cr + fi)sinht), 



Ly.jit) = [{j - <T){j + a-i + d)f^ y^-iAt), 



(A.15) 



The modes ^[/^'"'^ are explicitly 



(A.16) 



(A.17) 



^^""''^ = - j(j+7_i) (^°"^^) [{dt + {d-l)tanht)y,,j{t)\ ('^)l),Wyj(17rf) , (A.18) 

where is a normalization constant whose value is not important here. By inspection 
we see 

Vt(^oj = dtyoj{t)'^''^Yj{nd) 



[j{j + d-i)] 

—i 



1/2 Ar(*) 



A 



t ' 



(A.19) 



and 



N. 



it) 



^ju+d-1) A^^-^M 



(A.20) 



This proves that ^^^'■'^ oc V^c^q j. 



The modes A^f'^^ 

Although the modes A^ff'''^ should define an irreducible representation of the de Sitter 



group, the action of £5 on is not a linear combination of such modes. The issue 

;on 

term c£V^(/9o in (5_b (3.12) is non-trivial. 



is simply that £bA^i^'^^ also contains a pure gauge term which must be removed; i.e, the 
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We begin by examining the t and x components in detail: 



id-J) 



A 



X 



0, 



k{k + d-2) 
(i + l)(i + d-2) 



1/2 



cosh t 
sinx 



d-l) 



(A.21) 



First consider the boost action on a['^'''^] because a['^'''^ 
gauge contribution: 



it is easy to extract the pure 



6bA 



{d;j) 



smx 



. 



(A.22) 



cosh^ t ^ 

Next we use the a raising and lowering operators (A. 15) to note 

Ra=^iLa=oVo,j = -{coshtdt + {d - l)sinht)(cosht5j)(/?gj = j(j + d- l)v'oJ (A.23) 
(this can also be seen from (2.12)). We can now use this equation to invert (A.22), namely 

i(cosh tdt + {d-l) sinh t) ^1^4''''^ = + d- 1) VoJ ' 
from which we find 



(A.24) 



j{j + d - 1 



-V, 



(cosht 9t + {d-l) sinht) ^^4^^'^") 

cosh t ^ 



(A.25) 



We now have an explicit expression for the pure gauge contribution to the A^^'^^ modes. 
A.4 Tensor STSHs 

Finally, we turn to rank 2 STSHs on the sphere S"^. These harmonics ^'^^Tab satisfy 

{d)^ab id)j,^^ ^ id)jja id)j,^^ ^ id)^2 id)j.^^ ^ + d-l)+2] ("')r,fe. (A.26) 

The solutions depend on d angular momenta j which satisfy jd > Jd-i ^ • • • ^ i2 ^ 
\ji\ > 2. The {d + l){d — 2)/2 independent sets of solutions to (A.26) are have non- 
vanishing components tangent to the sub-spheres 5^, 5^, ... , S'^. The harmonics are thus 
labeled (d)T^f-^^^ where a, f3 run over the sub-spheres 3, . . . ,d. Explicit expressions for the 
harmonics are 

j.idd-J) ^ Nf,'\smxr'^''^yMx)^'~'%{^d-i) , 



rp{dd;j) 

rp{da-J) 
XX 

rp{da-J) 
X* 



rp{da;j) 



XX 



N 



(dd) 



jk 



k{k + d-2] 
0, a ^ d , 



[d^ + {d-2) cot x] ('^3^fc£(x) 



a ^ d , 



2m 



id) 



[k{k + d-2) -2] 



— (sinx)2 [d^ + {d-l) cot x] (")3^,fe(x) ("-i)l)(/"-^Vj"^'-)(17,_i), a^d, 



X« 



0, a, P ^ d , 



(A.27) 
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The normalization constants are 
(d) _ \{k-l){k + d-iy^'^ 



jk 



2jU + d-l) 



(dd) 



jk 



{d - 2){k - i)k{k + d- 2){k + d-iy^/^ 

{d-l)j{j + l){j + d-2){j + d-l) 



Tensor harmonics are normahzed with respect to the inner product 



(r«,r(2))^ := / di7,T«rM" = 6^'^^'\ 



(A.28) 
(A.29) 



B. Finite boost matrix elements 



In this appendix we compute matrix elements Bj^{X) = (^B{X)(p^j,if„^'j^^. This reduces 
to computing the integral 



Bjmk{X) = -2i dty^^it) I dx{smxy-^y.j{t')^''^yjk{x)'^''^ymk{x) 



(B.l) 



i=0 



where yaj{t) and ^'^^yjkix) a-re defined in terms of Legendre functions in (2.14) and (A. 2) 
respectively, and t' , x' a-i's defined in (2.21). The Legendre function may be defined in terms 
of the Gauss hypergeometric function via 

m/2 

2F1 



1 



r(i + m) 



z - 1 



z + 1 



-n , 1 + n ; 1 + m ; 



1 - z 



(B.2) 



After some standard manipulations of the hypergeometric function we may write 3^o-j(i) 
and ^'^^yjkix) as 



y{t)aj = Nd^jicoshty 



X Cf{smhty2Fi 
/=o,i 



^'-" + ^^' + " + ^ l+/;-sinh^t 
2 2 '2 



)y{x)jk = Mdjk{smxricosxy-'2Fi 



k — j k — j + 1 d 2 

, ^ ; /c + - ; - tan x 



where the constants N^aj, ^djk and Cf are 



Mdjk 



1 



r [j - j + a + d\ 



2i+(d-i)/2r [j + e^il) [ 
1 



1/2 



j + A; + d - 1 
j-k + l 



1/2 



2 J ■:J ' 2 

2 ' 2 



(B.3) 
(B.4) 

(B.5) 
(B.6) 
(B.7) 



Here r[. . . ] denotes products (and quotients) of gamma functions. Utilising these expres- 
sions and noting that 



dtyamit) 



t=0 



2i^J^^:Ndam^ 



m + 



d+l 



m~a m+a+d 
2 ' 2 



(B.5 
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the boost matrix elements (B.l) become 



where I is the integral 



771 + 



d+l 



m—a m+cr+d 
2 ' 2 



(B.9) 



^ = E C'/^''dx|(sinx)'^'^+'Xcosx)™-'(smxO'(cosx'F"'(cosht')^(sinht')^ 



/=o,i 



X 2F1 



X 2F1 



X 2F1 



k — mk — m + 1. d n 

' — 2 ; ^ + 2 ; - X 



fc — j A; — j + 1 d 2/ 



^'-^ + ^^' + " + ^ l + /;-sinh^t-' 
2 2 '2 



(B.IO) 



Let us now focus on /. First note that at t = the boosted coordinates satisfy (see 
(2.21)) 

sinht' = sinh A cos Xi sin x' cosh = sin x, (^ = 0) (B-H) 
cos x' cosh t' = cosh A cos X 5 tan^ x' = cosh^^ Atan^ Xi (* = 0)- (B.12) 



We may use these relations to eliminate the boosted coordinates from (B.IO): 

/ = (coshA)^'^^' ^ C7j(sinhA)-^ 
/=o,i 

^ [dx\ (sin x)'-'+''(cos ^y+^+f-^^ 



X 2F1 



X 2F1 



X 2F1 



k — m /e — m + 1 , d 



2 ' 2 

k-j k- j + 1 
2 ' 2 



; A; + - ; - tan x 



k + -] — cosh ^ A tan^ x 



j - + / i + o- + 1^ 1 , . . , 2 X 2 

— ^ — , — ^ — ; 2 + / ; ~ si'^^ '^os X 



(B.13) 



We can now expand the first two hypergeometric functions into their power series repre- 



sentations, i.e., 



2-t^l 



k — mk — m+1 d 9 

' 2 ■,k+- - - tan X 



[(m-fe)/2] 



^ ^^(-tan^x)". 



(B.14) 



a=0 
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where with condensed notation we take ipa to represent the coefficient of the power series 

/ k—m \ f fc— m+1 '\ 
^ 2 Ja\ 2 )a 



(fc + |)^r(l + a) 



(B.15) 



and (a)n := r(a + n)/T(a) is the Pochhammer symbol. With this / becomes 

[(m-fc)/2] [(i-fc)/2] C 

/=0,1 a=0 b=0 I 

d-l+2k+2a+2bf^^^^^y+ni+f-2k~2a-2b 

r j - + / j + 0- + d 1 2 X 2 ^ 

X 2-ri , ;?:+/; — smh A cos x 



X / dx(sinx)'^-^+^'^+^'^+^''(cosx) 
Jo 



(B.16) 



Because the integral is over the region to vr, the integral is non-vanishing only if the 
integrand has an even number of cos x- Thus of the two terms / = 0, 1, only one survives: 
/ = 0(1) if J + m is even (odd). 

The integral (B.16) can now be computed using the formula 



aJ'b 



a;b;zy 



c,d — c 
d 



a+iFb+i 



a , c ; 6 , d; z 



, (B.17) 



where a and b are lists of parameters of length A and B respectively. Thus / is 



Y ( j+m+f+d+l 



(coshA)^-*^(sinhA)-^ 

l(m-k)/2]l(j-k)/2] r 



a=0 6=0 



K + - + a + 6, ^ k - a-b 



(cosh A) 
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X 3-?^2 



j+cr+rf+Z j+m+/+l ,^ „ 
2 ' 2 ' 2 ft. u, u , 



1 I f j+m+d+l+f . 

2 ' 2 ' 



sinh^ A 



(B.18) 



Using this result we can now write the boost matrix elements as 

r [(m-fc)/2] [(i-fc)/2] 

B^mk{X)=A"^^k{cos\iXy-\ism\i\)n E E ^"™,fc(coshA)- 

I a=0 6=0 



26 



X 3-?^2 



2 ' 2 ' 2 ft. u y , 



1 I f j+m+rf+l+/ . 

2 ^ J 1 2 ' 



- sinh^ A 



(B.19) 
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where explicitly 



1 



22k+d-2-fY ( j+m+f+d+1 




1/2 



(B.20) 




— 1 \ 

— ) {k- m)2a {k - j)2b r 



(B.21) 
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